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Abstract 



A signed measure analysis of two-dimensional intermittent magnetohy- 
drodynamic turbulence is presented. This kind of analysis is performed to 
characterize the scaling behavior of the sign-oscillating flow structures, and 
their geometrical properties. In particular, it is observed that cancellations 
between positive and negative contributions of the field inside structures, are 
inhibited for scales smaller than the Taylor microscale, and stop near the dissi- 
pative scale. Moreover, from a simple geometrical argument, the relationship 
between the cancellation exponent and the typical fractal dimension of the 
structures in the flow is obtained. 
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I. INTRODUCTION 



Self-similarity is a signature of complex flows with strong nonlinearities. Huge amount of 
efforts have been deployed in determining scaling laws of energy spectra and, more generally 
speaking, of structure functions (see e.g. Ref. 1 and 2), both in geophysical and astrophys- 
ical flows, in the laboratory and using numerical simulations as well. Non-linear behavior 
of the scaling exponents is generally interpreted as due to the presence of strong localized 
structures both in space and time. Indeed, structures have been observed in turbulent flows 
for a long time: sheets, spirals and filaments of vorticity both in the incompressible 
and compressible case 0, in the latter with shocks as well. In magnetohydrodynamics 
(MHD), current and vorticity sheets seem to prevail in three dimensions [PHT^ whereas 
in two-dimensional space, current and vorticity filaments are the main identifiable struc- 



tures ||1^H16|. Such features are locally complex, with rapid fluctuations corresponding to 
flux cancellations, in three dimensions [|10[ as well as in two dimensions, and as first strik- 



ingly observed in chaotic dynamos using mappings |]T^; in such cases, very high Reynolds 
number flows can be modeled and cancellations of flux occur at all scales, leading to a very 
complex magnetic field. 

We study in this paper the geometrical properties of the vorticity and current fields (as 
well as their combination through the Elsasser variables defined below). The data stems 
from a direct numerical simulation in two space dimensions and our study focuses on the 
computation of the cancellation exponent introduced by Ott et al. [|l^. These exponents 
allow for a simple characterization of the flows and can be linked to the fractal dimension 
of the typical structures. 

The present paper is organized as follows. In Section we recall the definitions of the 
signed measure and cancellation exponent. In Section |T|, we describe the data stemming 
from a numerical simulation, and the main features of the fields are characterized by the 



typical scales of the flow. Section |^ describes the results obtained from the analysis of the 
cancellation properties of our flow. In Section |V|, we introduce a simple model, based on 
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geometrical arguments, which connect the cancellation exponent to the fractal dimension 
of the structures. Our results are discussed and the main conclusions of the paper are 
re-emphasized in Section |VI[ 



II. THE SIGNED MEASURE AND THE CANCELLATION EXPONENT 

Intermittency in turbulent fields is related to the presence of structures at different scales. 
The proximity of such structures may lead to rapid changes of sign in the fields derivatives. 
In order to describe the scaling properties of sign oscillations, Ott et al. [|l^ introduced the 
concept of sign-singular measures. In analogy to probability measures (positive definite), 
the signed measure of a zero-mean scalar field /(r), defined on a set Q{L) of size L, can be 
introduced in the d- dimensional case. Let {Qii})} C Q{L) be a hierarchy of disjoint subsets 
of size / covering Q{L). Then, for each scale I and for each set of boxes Qi{l), the signed 
measure is defined as 

^'^ ^ hiL, dr \fir)\ ■ 
The choice of the normalization factor (the denominator in Eq. (|I])) will be justified be- 
low. The signed measure thus can be interpreted as the difference between two probability 
measures, one for the positive part and the other for the negative part of the field JTSf. 

As the scale of the subset Qi{l) increases, cancellations between opposite sign struc- 
tures, and then between the two probability measures, come into play. This feature can be 
characterized by investigating the scaling properties of the partition function : 

x(0=E 1/^^(01 (2) 

Q^{l) 

where the sum is extended to all disjoint subset Qi{l). It is clear that if /ij(/) is a probability 
measure, as well as if the field is defined in sign, the choice of the normalization factor in (|l]) 
leads to x(0 = 1- other hand, if cancellations are present, we can expect that their 

effect is less and less important as the scale / decreases, so that x(0 increases. In this latter 
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case, we can investigate the scaling behaviors of the weakness of cancellations by a scaling 
exponent k, defined through ~ where k is called the cancellation exponent. This 
exponent represents a quantitative measure of the cancellation efficiency. Furthermore, if 
this scaling law exists, in analogy with the singularities of probability measures, the signed 
measure is called sign-singular. If the field is smooth, then trivially k = 0. It can be shown 
that for homogeneous fields with discontinuities, k = 1, whereas for stochastic processes 
K, = d/2 [jl9|, d being the dimensionality of the space. More generally, as will be shown 
in section [V^ , if the field f{r) is homogeneous with a Holder scaling exponent h, that 



is if (||A/(Z)||) = (||/(r + Z) - /(r)||) ~ V\ then k = -h [|T|,|g]. This relation between 
the cancellation exponent and scaling exponents can be generalized to higher order partition 
functions and structure functions S'p(/) = (AII/(/)p) ~ I'^p (where All/(Z) = [f{r+l) — f{r)]-l 
defines the longitudinal increments of / on scale /) to determine the generalized dimensions 
of the set where the singularities of the field / live |[T9|j2T[|. 



Besides the work of Ott et al. on fast magnetic dynamos, several fiows in the context 
of plasma turbulence have been examined using such concepts, in order to characterize the 
oscillating features of turbulent fields. For example, the magnetic helicity in the solar wind 
has been shown to be sign-singular , and solar photospheric velocity patterns as well. 
Finally, the variations of the cancellation exponent of the current helicity between the pre- 
and post-fiare periods for major H-a solar fiares, have been used to evidence variations in 



the heliospheric magnetic field structures of active regions |24 



III. THE NUMERICAL SIMULATION AND THE TYPICAL SCALES OF THE 

FLOW 

The data we analyze with the technique described in the previous section stems from a 
numerical simulation of the two-dimensional {d = 2) magnetohydrodynamics (MHD) equa- 
tions written below in terms of the Elsasser variables = v ±b 

[dt + ■ V)z± = -VP* + z/^V^z± + z/^V^zT + F± (3) 
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where is the total pressure, 2z/^ = u ± rj, v being the viscosity and 77 the magnetic 
diffusivity, while are the forcing terms. The incompressibility conditions V ■ = 
complete the equations. 

The MHD equations have nonlinear scaling properties so that Kolmogorov-like argu- 
ments ||2^ can be performed to evaluate dimensionally the scaling properties of structure 



functions, but taking into account the specificity of MHD flows. The less efficient interaction 
rate between the structures may well lead to scaling laws that differ from the classical Kol- 
mogorov case, as first derived in the so-called Iroshnikov-Kraichnan phenomenology (IK) ||2^ 
for globally isotropic flows {i.e. in the absence of a strong uniform magnetic field Bq at zero 
frequency), a phenomenology which takes into account the presence of Alfvenic fluctuations 
in the plasma. Small-scale intermittency is present in MHD as in the fluid case, and is 
related to the presence of coherent structures which locally in space obviously break the 
isotropy assumption. However, in the Kolmogorov spirit, isotropy is recovered on average. 
The effects of these structures on the phenomenology for MHD flows have been pointed out 
and discussed in several ways (see, e.g. Ref. 10 and 27). 

In the present paper, we want to analyze the oscillating character of such structures, and 
their geometrical properties. Equations (|^) are solved with a pseudo-spectral method on a 
grid of 1024^ points, with 27r-periodic boundary conditions. The forcing terms maintain 
constant the Fourier modes k with \k\ = 1. The magnetic Prandtl number = v/rj is unity, 
so that z/+ = z/ = ?7 = 8x 10~^. The correlations between the velocity and the magnetic 
fields leads initially to PC = 2('y-b)/(|'i;p + |6|2) ^ Q%. The time -averaged Reynolds number 
is Re ~ 1600 and the integral scale is Iq/ L = 0.25 ± 0.02, where L = 27r is the size of the 
computational box. Further details about the simulation can be found in Ref. 28 and 29. 
In the present work, we analyze ten samples separated by approximately 16 eddy turnover 
times once all transients have died out and a statistically steady state has been reached. 

A snapshot of the current is shown in Figure |1|. The top picture shows the local ac- 
cumulation of sheets, probably through folding, whereas the bottom picture clearly reveals 
the oscillating character of the flow structures. The analysis of the probability distribu- 



tion functions (PDFs) of the field fluctuations and of the scaling exponents of the 
structure functions I^Hf has revealed that both the Elsasser fields, as well as the magnetic 
field and the velocity, are strongly intermittent, as already visualized in earlier computa- 
tions This intermittency is already quite evident when one computes the flatness 
factor F{s) = ([AII/(s)]'')/([AII/(s)]^)^ that goes from the gaussian value 3 at large scales 
s for all fields, to values F ~ 40 for / = i?, F ~ 80 for / = b, F ~ 65 for / = at the 
smallest resolved scale. This dramatic change in flatness reflects the change in shape of the 
PDFs through the scales, as quantified in Ref. 29. 

Following Politano and Pouquet |3^, it can be shown that the third order correlators 
defined below obey, in the inertial range, the following scaling law in dimension d 

y3±(0 = (A"^^(0l|A^±(0lP) = -^e±/ (4) 

where Az'^^l) = z^{r + Z) — z^{r), is the vector increment of the field z^, and A" is its 
(scalar) longitudinal projection; thus, the flux functions Y^{1) involve all components of the 
physical fields. The coefficient are the mean energy transfer rates of the z^ variables. We 
recall that in the framework of the Kolmogorov phenomenology, the mean energy transfer 
rate is assumed to be equal to the mean energy injection rate, as well as to the mean energy 
dissipation rate. In order to estimate e^, the data are fitted for each time in the ranges 
where the relation is verified. Averaging in time in order to decrease statistical errors, 
we evaluate ~ 550 and ~ 440. The difference between the two ± rates is linked 
to the correlation between the velocity and the magnetic field, and should increase with 
increasing pc and disappear only in the limit of vanishing pc- Figure || displays the Y^{1) 
third order correlators. It appears that the exact scaling law (^) - obtained in the limit of 
very large Reynolds numbers, otherwise correction terms not written here would arise due to 
dissipation - is better verified in the case of 1^3^ than for Yg", where the scaling range coincide 
with the inertial range as determined in Ref. 28 (see below). Note that, because the fluxes 
are odd-order in the fields, cancellations occur that render scaling laws difficult to unravel. 
For that reason, it is customary to estimate the inertial range by computing the modified 



third order correlators, taking absolute values to obtain Lf(l) = (|Al'z'^(/)| ||Az^(/)||^) ~ /. 
This leads to an inertial range extending from 0.01 to 0.1 'ml/ L for both correlators p8[. The 
difference in scaling quality between the plus and minus third order correlators, or higher 
order moments of the 2;+ and z~ fields, has already been observed, even when using absolute 
moments to reduce statistical errors This asymmetry between and z~ could 

be due to a significant role of the velocity-magnetic field correlations on the dynamics of 
MHD flows, even when ~ 0. Indeed, the local v ■ b, not positive definite, is known to 
develop strong fluctuations (see, e.g. Ref. 31), precluding identical behavior for fields. 

We can use the values of computed from the law (^) to estimate - in the framework 
of the IK phenomenology - the typical scales of the flow, i.e. the dissipative scales 

= {u'^BQ/e'^)^^^ and the Taylor microscales = (/luE^ /e^)^^'^ . Here Bq represents the 
r.m.s. large-scale magnetic field, and are the energies associated with the Elsasser fields, 
which are the ideal invariants of the flow, together with the (a^) correlation, where b = V xa 
with a the magnetic potential. The computed scales we obtain are given in Table |, and 
indicated in the orizontal axis of Figure ^ as well. In MHD turbulence, it is not possible to 
unambiguously define directly from the dissipation rates, the Taylor and the dissipative scales 
for the magnetic and velocity fields independently. However, using several assumptions, 
namely (i) equipartition, i.e. {v'^)/2 ~ (b^) /2 for the kinetic and magnetic energies, (ii) weak 
total correlation, (1; • 6) ~ 0, and (iii) unit magnetic Prandtl number, ~ 1, one can then 
define a dissipation scale Id = (z/^-Bo/e)^^^ and a Taylor microscale A = (2z/i?/e)^^^, where E 
stands for the total energy, and e for its mean transfer and dissipation rates. Doing so, and 
noting that e = (e+ + e-)/2, we obtain la/L = (2.2 ±0.4) x 10~^ and X/L = 1.5 ±0.4 x lO'^. 
These scales are very close to the ones computed for the Elsasser fields. We are thus left to 
use and A^ as the typical scales of the flow to refer to. 
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IV. SCALING PROPERTIES OF THE SIGNED MEASURE IN THE 

NUMERICAL DATA 



We now analyze the scaling properties of the signed measure for the numerical data 
described above. As in the previous section for the evaluation of the energy transfer rates, 
we do so separately for the different snapshots of the simulation, and we report only the 
computed averages, while the error bars give an estimate of the dispersion of such values. 

We first define the signed scalar densities we need from the physical variables b 
and z^. For two-dimensional MHD flows, it is convenient to use the rotational of these 
variables. Indeed, using the Stokes theorem, the circulation of the magnetic field along 
any closed contour, is equal to the normal flux of the current crossing the corresponding 
surface Qi{l) §Ciii)^i''^) ' ~ /q,(/) '^('") ■ n da . The sign of the current density flux is 
determined by the clockwise or counterclockwise circulation of the magnetic field. In a two- 
dimensional geometry, the magnetic field lies in the (x, y) plane, namely b{x, y) = {b^, by, 0), 
and the current density J{x,y) = V x b{x,y) is along the z-direction, perpendicular to 
the plane. The normal current flux density is then simply J ■ hda = Jz{x,y) dx dy, so 
that Jz{x, y) = J{x, y) can be used as our scalar density. The same argument holds for the 
vorticity uj and for uj^ = uj :k j, the rotationals of the Elsasser fields. We show in Figure |^ 
the coarse-grained current density flux for four different values of the box size /, at one 
particular time (t = 7.3 ~ 45 t^vl, where t^l is the eddy turnover time once the statistically 
steady state is reached). As can be seen, strong oscillating signed structures appear at all 
scales, and cancellation effects clearly increase with the box size. In Figure ^ we present the 
function X^i})^ computed in the case of the current J{x,y) at each time and averaged over 
the ten temporal fields we analyze. A clear scaling range is visible, so that a cancellation 
exponent = 0.43 can be obtained by performing a least-square fit of the type 

X(0 = C \ (5) 

We now try to relate the scales which characterize the cancellations to the typical scales of 
the flow itself, namely the Taylor microscales and the dissipative ones l^. The scaling 
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range is determined by looking at the local slope of (not shown), and extends from 

= 0.015L to I'l = 0.12L. As can be observed, the range in which the scaling (^) holds 
is mostly embedded within the inertial range of the flow, given in section For scales 
larger than /„, the partition function still decays as a power-law, but with a larger exponent, 
approaching the typical behavior expected for a completely uncorrelated point field, with 
K = d/2 = 1, but the poor statistics at such large scales prevent us from discussing that 
behavior any further. For scales smaller than the partition function slowly saturates at 
the asymptotic value x(0 = 1- The scale Is, at which the saturation can be considered 
achieved, is computed as the intersection scale between the power-law @ and the limit 
value x(0 = 1) giving Ig = 0.005L for the current. It is worth pointing out that, due to the 
choice of normalization in the signed measure, for any dataset with a finite resolution /min, 
the partition function will trivially go to unity at that scale. The fact that, in our case, 
saturation of x(0 persists for scales larger than /min, demonstrates that the field is well 
resolved and thus regular at the resolution of the numerical simulation. The departure 
from a constant behavior of x(0 (with /t = 0), occurs only at some yet larger scale, which 
can be identified as the smallest characteristic scale of the structures which develop in the 
flow. Typically, we find that the saturation scale Is is of order of 2 3 times the relevant 
Kolmogorov scale Id, showing that the extent of structures in their smallest dimension is 
controlled by dissipative effects. On the other hand, power-law scaling sets in at the scale 
which is found to be of the same order as the Taylor microscale A. This point can be 
understood by a simple geometrical argument developed in the next section. All these results 
obtained for the current are confirmed for the other fields u, uj~^ and uj~ . The values of the 
cancellation exponents and those of the typical scales related to the partition function for 
all fields are reported in Table |. It can be observed that the rotationals of the Elsasser 
fields both present the same behavior, while different values of the parameters are found for 
the vorticity and the current. Note also that the range of fit for the vorticity is shorter than 
those for the others fields (not shown). 

The higher value of k'^ reveals that the presence of structures inhibit the cancellations 
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in a less efficient way for the vorticity than for the other fields. Conversely, the current 
structures inhibit cancellations more efficiently than all other fields. 



V. CANCELLATION EXPONENTS AND DIMENSIONS OF STRUCTURES 

It is clear that cancellation exponents characterize in some way the geometry of struc- 
tures, as exemplified by the fact that a smooth (continuous) field will give a zero cancellation 
exponent, while a field made of uncorrelated point-like objects will have a cancellation ex- 
ponent K — d/2. In order to quantify the transition between these two limits, we will use a 
simple geometrical model that we now introduce. 

The relevant Taylor scale A gives us a mean scale over which a field is correlated. How- 
ever, it does not tell us the geometry of the structures present at that scale. These could be 
correlated in some directions for scales much larger than A, and/or already uncorrelated in 
some other directions. Let us now assume that the field is smooth (correlated) in D dimen- 
sions with a cutoff scale of A, and uncorrelated in the other d — D dimensions. Of course, 
if the field becomes completely regular below some scale Ig, the corresponding dimension 
at that scale will he D = d. On the other hand, if structures have some largest extent /„, 
the field will consist of isolated, and thus nearly uncorrelated, objects, showing an apparent 
dimension D = In the intermediate region, over which the topology of the structures will 
have an effect, the partition function of the current J for example, or of any component of 
a field in dimension d, can be written as 



where homogeneity is used to replace the sum over all subsets Qi{l) by {L/lY times the 
integral over a generic box Q{1) of size I. Moreover, we have approximated the field abso- 
lute value integral by its typical value L'^J^is- Indeed, because of the absolute values, no 
cancellations can occur in the denominator. Now, the integral on the subset Q{1) can be 




(6) 
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done by integrating over regular domains of size A*^ and splitting the number of contribu- 
tions between the correlated dimensions and the uncorrelated ones. The smooth dimensions 
will give a contribution proportional to their area {l/X)^. The uncorrelated dimensions will 
behave like the integral of an uncorrelated field, giving a value proportional to the square 
root of their area, that is (1/ XY^^^^^"^ . Collecting all contributions, we see that will 
behave as 




The cancellation exponent can thus be interpreted as half the codimension of the flow struc- 
tures, K = {d — D)/2. Note that this interpretation remains valid at small scales where k = 
when D = d and at scales larger or similar to the integral scale where k = d/2 when D = 0. 

Specializing to our case with d = 2, the value of = 0.5 that we found indicates that 
structures are similar to filaments = 1.0 ± 0.12 for the o;^ fields. This result for the 
Elsasser variables rotationals extends the previous observation of sheet-like structures for 
the current in 3D MHD [jl0|,|12|,|13[ , having a corresponding signature of filaments in 2D cuts. 

The current displays a slightly smaller cancellation exponent, that leads to a dimension 
of 1.14 ±0.12, which can be interpreted by observing that the current structures are slightly 
fatter than filaments and show some thickness, likely due to reconnection taking place within 
current sheets. 

The vorticity field has a significantly larger cancellation exponent, that would give a 
dimension for structures D'^ = 0.62 ± 0.24, smaller than the one of a filament structure. In 
fact, a close inspection in physical space shows that vorticity structures, although globally 
thin and elongated, are more complex than those of the current field (see plots at t = 
6.3 ~ 26 TiVL, and at t = 6.93 ~ 39 t^l in Figure |]). For example, as already known ||3T1 , |32| , 
quadrupole vorticity structures are associated with quasi-linear current sheets in the simplest 
reconnection configuration. This could explain the smaller dimension found for the vorticity 
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structures characterizing what can be seen as a "disrupted filament" . 



VI. CONCLUSION AND DISCUSSION 



In order to characterize the oscillating behavior of the dynamical structures present in 
conductive turbulent flows, we perform a signed measure analysis of data obtained from 
direct numerical integration of the two-dimensional MHD equations. We obtain that the 
measure is sign singular, with a clear scaling range for the cancellations between the struc- 
tures, either for the current and the vorticity field, or for the rotational of the Elsasser 
variables. The cancellations are inhibited for scales smaller than the Taylor scales while 
they stop near the dissipative scales. Note that the correspondence - as displayed in Table | 
- between the typical scales for cancellation and those of the turbulent fluid are obtained 
in the framework of MHD flows presenting a weak correlation between the velocity and the 
magnetic fields and for a magnetic Prandtl number of unity. 

Moreover, by a simple geometrical argument, we link the cancellation exponent evaluated 
here to the fractal dimension of the structures in the flow. This gives us informations about 
the topology of the structures. This interpretation can in fact also be extended. Indeed, 
the partition function is only a first-order moment, and is thus unable to distinguish, for 
example, between a regular behavior over some fractal set, and a self-similar behavior over 
the whole space, as well as intermediate situations. The distinction could only be made 
by looking at higher order moments [0. This point is left for future study. But we can 
already generalize the geometrical argument given in the previous section to the case where 
the field is not locally smooth, by introducing a single local scaling exponent h such that 
(specializing to the magnetic field) ||AJ(s)|| ~ s^~^ on the structures of dimension D. Note 
that, because of the presence of structures, h is not the exponent of the first-order structure 
function of J, the latter containing contributions both of structures and background. In this 
picture, the last integral in expression (|^) becomes 
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D+h~l / , X ^ 



^ <-'rms 



2 



giving a cancellation exponent k = (rf — D)/2 + (1 — /i). In general, as we already pointed 
out, it is impossible to separate the contributions of the two terms, unless one makes further 
hypotheses. 

Let us consider the particular case of a space filling object (D = d) which would not be 
smooth at small scales, that is with a Holder self-similarity exponent h smaller than one. 
The cancellation exponent of its rotational would then be simply k, = 1 — h. Note that if 
the object is smooth, with h = 1, this is consistent with k = obtained from (J^). On the 
other hand, if the field values are completely independent, we find n = 1, that can be easily 



checked by direct numerical simulations |23|. Using this interpretation, in the case of the 
velocity rotational that has a cancellation exponent k = 0.69 ± 0.12, we are led to infer an 
Holder exponent for the velocity h = 0.31 ± 0.12, close to the Kolmogorov value 1/3 for the 
velocity field of a neutral turbulent flow. This argument is clearly incomplete because of the 
presence of structures in the fields, which are even more conspicuous for the rotationals of 
the magnetic and Elsasser fields. Higher resolution simulations allowing the computation of 
the full spectrum of higher order cancellation exponents will be necessary to disentangle the 
contributions of geometry and differentiability to the intermittency of MHD turbulence. 
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TABLES 

TABLE I. Estimated values for the typical scales of the flow (see text for definitions), the 
cancellation exponents and the fractal dimensions of the striictures for the four fields, averaged on 
ten temporal snapshots; la is the dissipation scale, A the Taylor scale, Z* denotes the end of the 
inertial range and Ig is the saturation scale. 
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FIGURES 

FIG. 1. A snapshot of the current J at time t = 7.3 in grey levels (top), and a pseudo-3D 

perspective view of the same field (bottom) showing the existence of structures of different signs 
at all scales. 
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FIG. 2. The time averaged third order correlators Y^{1) (continuous hne) and Y^{1) (dashed 
hne) in a log-log plot. The straight line is a simple linear law, and is displayed as reference. The 
dissipation scale Zj", the Taylor microscale A"*" and the integral scale Iq are indicated, as well as the 
bounds of the inertial range. 
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FIG. 3. The coarse-grained signed measure of the current J at time t = 7.3 for four different 
box sizes, namely l/L = 0.001, l/L = 0.016, l/L = 0.059, l/L = 0.12, from top to bottom. Colors 

range from cyan for negative J values to yellow for positive ones, going through blue and brown. 
Cancellations at large scales are responsible for the decrease in magnitude of the measure. 



19 



10' 









l/L 


I'/L 





10' 10" 10' 10° 



UL 

FIG. 4. For the current J, the log-log plot of x(0 (averaged over ten snapshots) versus the 
subsets size l/L. The solid line represents the best fit with the power-law (^). The constant value 
Xs(/) = 1 is plotted (the horizontal dashed line), and the saturation scale Is/L = 0.005 is indicated 
by the vertical dashed line on the left. The vertical dotted lines indicate the range of the fit, which 
lies between U/L = 0.015 and h/L = 0.12. 
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FIG. 5. Examples of contour lines of vorticity (left panel) and current (right panel) structures, 
at the time t = 6.3 (top) and t = 6.93 (bottom). The axes indicate the location of the snapshots 
within the 27r-periodic box. 
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